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Four-neutrino mixing solutions of the atmospheric neutrino anomaly 

A. Marrone'^ 

'^Dipartimento di Fisica and Sezione INFN di Bari 
Via Amendola 173, 1-70126 Bari, Italy 

Solutions to the atmospheric neutrino anomaly which smoothly interpolate between i/^ — > Vr and 
oscillations are studied. It is shown that, although the Super-Kamiokande data disfavor the pure v^^ 
channel, one cannot exclude sizable amplitude for the — > Vs channel in addition to Vr oscillations. 
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1. Introduction 

Four-neutrino (4z/) [|]-|| models including a hy- 
pothetical sterile state (r^s) can accommodate the 
three sources of evidence for v flavor oscillations 
coming from atmospheric, solar, and accelerator 
neutrino experiments. In particular, \v spectra 
with mass eigenstates organized in two doublets 
(2-1-2 models) seem to be favored by world neu- 
trino data 1^. Although 2-1-2 models are often as- 
sumed to imply the v^, — > Vr or i/^ — > Vs channel 
for solar and atmospheric oscillations, one might 
have mixed (active-|-sterile) flavor transitions of 
the kind S 



Vj^ (atmospheric) , 
(solar) , 



(1) 
(2) 



where the states v±^ as discussed in detail in 
represent linear (orthogonal) combinations of Vr 
and Vs through a mixing angle 5. The oscilla- 
tion modes (|l|,^) represent generalizations of both 
modes v^J^ — > Vr and — > i/^, to which they re- 
duce for sinf = and 1. For generic values of 
sin^, the final states in atmospheric and solar 
j/g flavor transitions are linear combinations of 
and i/r, and the coefficients have to be constrained 
by experiments. In this work, four-neutrino oscil- 
lations in the context of 2-1-2 spectra are studied, 
for unconstrained values of sin^ = (i/_|_|z/s). It 
is shown that the state (into which oscil- 
lates) can have a sizable component. As shown 
in ^ , these results can be reconcilied to the Av so- 
lutions to the solar neutrino problem [Qjsj, which 
are compatible with a large Vs component of v- . 



2. Graphical representations 

The mixing parameter spaces for atmospheric 
and solar neutrinos can be represented in tri- 
angular plots, embedding unitarity relations of 
the kind Ul + Ui + Ui = 1, Q holding for the 
four columns of the mixing matrix U, under suit- 
able approximations It can be shown that 
is sufficient to implement just one unitarity rela- 
tion in one triangle plot for the two mass dou- 
blets, the first (1/3, 1^4) being the atmospheric 
doublet, the second (1^1, 1^2) the solar doublet, in 
our framework. Under the semplifying hypothe- 
ses discussed in Q, the does not take part to 
the oscillations of atmospheric neutrinos, so that 
the mixing can be graphically represented in a tri- 
angle plot, where the upper, lower left and lower 
right corner are identified with f^, z/j, and v^-, re- 
spectively, and the actual mixing is determined 
by the position of the 1/4 mass eigenstate in the 
triangle. The heights projected from a generic 
point i>4 inside the triangle onto the lower, right, 
and left side represent the elements U^^, and 
respectively. These three matrix elements 
can be expressed as functions of two mixing an- 
gles (-0 ,^). The square mass difference be- 
tween the atmospheric doublet (^3, 1/4) completes 
the parametrization. When 1/4 coincides with one 
of the corners no oscillation occur. Generic points 
inside the triangle describe mixed (active-|-sterile) 
atmospheric neutrino oscillations, smoothly inter- 
polating from pure v^ ^ s to pure ^ Vt-. 



^Triangle plots have been already introduced 
cussed in detail in the context of Si' mixing, see 
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Figure 1. Super-Kamiokande data on zenith 
distributions of lepton events induced by atmo- 
spheric neutrinos. 



3. Results of the atmospheric v analysis 

The most recent 70.5 kTy Super-Kamiokande 
data used in our analysis ^ are shown in Fig- 
ure 1. The statistical analysis is similar to 
that of jl^. Figures 2 and 3 show the results 
of our analysis in the atmospheric triangle plot. 
The first three columns of triangles refer to the 
separate fits to sub-GeV electrons and muons 
(10-1-10 bins), multi-GeV electrons and muons 
(10-1-10 bins), and upward stopping and through- 
going muons (5-1-10 bins), while the fourth col- 
umn refers to the total SK data sample (55 bins). 
For each column, we find the minimum from 
the fit to the corresponding data sample, and 
then present sections of the allowed volume at 
fixed values of for Ax^ = 6.25 (90% C.L., 
solid lines) and Ax^ = 11.36 (99% C.L., dot- 
ted lines). The low-energy SG data are basically 
insensitive to s^, since they are consistent both 
with j/^ Vt oscillations (left side) and with 
— > Vg oscillations (right side), as well as with 
any intermediate combination of the two oscilla- 
tion channels. High-energy upgoing muon data 
are instead much more sensitive to the compo- 
nent through matter effects, which increase both 
with energy and with s|, and tend to suppress 
the muon deficit. Large matter effects appear 
to be disfavored by the upgoing muon data, the 
rightmost part of the triangle being excluded at 







Figure 2. Results of the global fit to atmospheric 
neutrino data in the triangular plots, for de- 
cresing vn} values. 



90% C.L. The multi-GeV data cover an interme- 
diate energy range and are not as constraining as 
the upgoing muons. However, they show a ten- 
dency to disfavor a large sterile component, that 
is strengthened in the global combination of data 
(fourth column), leading to the upper bound 

s\ < 0.67 (90% C.L.) , (3) 

which disfavors pure or quasi-pure — > oscil- 
lations. From Figs. 2 and 3 we derive that pure 
i^fj. — > J^s oscillations are excluded at > 99% C.L., 
in agreement whith SK results 0. However, the 
bound (H) still allows intermediate cases of oscil- 
lations with a sizable I's component. For instance, 
one cannot exclude that the sterile neutrino chan- 
nel may have the same amplitude as the active 
one (s| — 0.5), or may even be dominant. Fur- 
thermore, for in its upper range, a subdomi- 
nant Vs component actually helps in the global fit 
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Figure 3. As in Fig. 2, for fower values. 



to SK data, leading to a 90% C.L. allowed region 
which does not touch the left side of the triangle 
(pure — > Vt channel). 

Figures 2 and 3 also provide bounds on the 
other mixing parameter s^, which governs the 
relative amount of and on the atmospheric 
neutrino states and and thus the overall 
amplitude of — > v+ oscillations. The following 
limit can be derived 



0.51 ±0.17 (90% C.L.) , 



(4) 



which favors — *■ oscillations with nearly 
maximal amplitude, as expected from the ob- 
servation of nearly maximal average suppression 
50%) of upgoing MG muons. Therefore, on 
the basis of present SK data on the zenith distri- 
butions of leptons induced by atmospheric neu- 
trinos, pure — > Vs oscillations are strongly 
disfavored, but one cannot exclude mixed active- 



4. Summary and conclusions 

In the context of 2-f 2 neutrino models, the SK 
atmospheric data have been studied, assuming 
the coexistence of both — > Vr and ~> Vs 
oscillations, with a smooth interpolation between 
such two subcases. It has been shown that, al- 
though the data disfavor oscillations in the pure 
— > Us channel, one cannot exclude their pres- 
ence in addition to oscillations. High en- 
ergy muon data appear to be crucial to determine 
the relative amplitude of the active and sterile os- 
cillation channels for atmospheric neutrinos. 
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